American Journal of Biological Anthropology

Extremely elongate necks have convergently evolved in several amniote lineages, including both aquatic and terrestrial forms (Fig. 1). The development of such a feature brings with it advantages in obtaining food items, but also biomechanical challenges, such as flexibility, stability, lift, and inertia. In Tanystropheus, a particularly long-necked Triassic archosauromorph, the neck is composed of only 13, mostly extraordinarily elongated and slender cervical vertebrae and accompanying rod-like, overlapping ribs, making it arguably the most extreme example of neck elongation in tetrapod evolution (Fig. 1;1-6). Understanding the function of this remarkable neck provides insights into the limits of neck elongation in amniotes and the evolution of morphological novelties in Triassic reptiles. Here we present the first quantitative biomechanical analysis of the Tanystropheus neck using a digital model based on three-dimensionally preserved bones. We assessed its range of motion (ROM) and performed finite element analysis (FEA) on the individual cervical ribs and the neck model in different configurations. Our results indicate that the neck of Tanystropheus was not extremely stiff, as previously postulated, and the ribs likely did not impair its movements. They transferred tensile forces towards the base of the neck, similar to what hypothesized for sauropods7. This study elucidates the bauplan of an extremely specialized animal and brings us closer to understanding the patterns of achieving neck elongation in vertebrates.

Birds evolved large, cognitively capable forebrains independently of mammals, yet comparative analyses of avian brain organization have been constrained by the lack of standardized resources capable of resolving internal parcellation and long-range connectivity across species. Here, we present a comparative MRI resource spanning 16 avian species representing major clades and diverse ecological niches. We analyzed high-resolution T2-weighted and diffusion-weighted datasets suitable for direct interspecific comparison. T2-weighted morphometry revealed pronounced region-specific variation in internal brain architecture, including lineage-dependent differences in the relative prominence of major brain divisions and commissural structures, supporting a pattern of mosaic diversification rather than uniform scaling. To validate MRI-derived anatomical boundaries, we compared MRI parcellations with complementary histological analyses in three representative taxa (the large-billed crow, gentoo penguin, and mandarin duck), demonstrating close correspondence between MRI-defined borders and cytoarchitectonic transitions identified by Nissl staining, as well as major myelinated compartments visualized by Luxol Fast Blue staining. Moreover, diffusion MRI tractography and fractional anisotropy (FA) mapping further revealed both conserved and species-specific features of large-scale brain organization. Seed-based tractography of the optic lobe, dorsal cortex, cerebellum, and anterior cortex in chick, gentoo penguin, and large-billed crow revealed conserved within-compartment trajectory patterns alongside marked region-specific interspecific differences, particularly in optic-lobe-associated long-range trajectories. Whole-brain FA maps revealed complementary variation in regional microstructural organization across taxa. Together, this comparative MRI framework provides a cross-validated foundation for linking internal brain anatomy and long-range connectivity to ecological and evolutionary diversification in birds, with broader applications to comparative neuroanatomy across amniotes.

Sharks exhibit extraordinary morphological diversity across a wide range of ecological niches, yet large-scale, high-resolution digital datasets of their internal anatomy remain limited. Here we present an open-access 3D shark anatomical repository derived from published X-ray computed tomography (CT) data, featuring manually segmented and systematically annotated models of the chondrocranium, visceral arches, axial skeleton, musculature, and viscera in standard STL format. The dataset comprises 117 individuals, representing 72 species across 25 families and all nine extant shark orders, with 115 full-body reconstructions and two head-only models. This open-access dataset offers a comprehensive resource for comparative anatomy, biomechanical simulations, evolutionary developmental biology and biomimetics research of extant sharks.

Adrenarche, the pre-pubertal rise in adrenal androgens, particularly dehydroepiandrosterone (DHEA) and its sulfated form DHEAS, is a critical driver of middle childhood cognitive and social development in modern humans. Compared to other apes, modern human adrenarche is more prolonged, with higher DHEAS levels. Whether this uniquely prolonged human adrenarche is a derived trait of Homo sapiens or has deeper hominin roots remains unresolved. Here, we examine the Neanderthal genetic variation in five key DHEAS biosynthesis genes (HSD3B2, CYP17A1, POR, CYB5A, SULT2A1). We also examine archaic introgression in these genes by comparing high-coverage Neanderthal genomes with globally diverse modern human sequences from the 1000 Genomes Project. We identify 29 Neanderthal-derived single nucleotide variants across these genes. Key steroidogenic genes in the biosynthesis pathway show no evidence of introgression, consistent with selection on pleiotropic regulators of steroidogenesis. In contrast, accessory genes carried introgressed Neanderthal haplotypes at moderate frequencies in non-African human populations, indicating Neanderthal variants are compatible with the human DHEAS synthesis pathway. All Neanderthal-specific variants were in non-coding regions, with three variants associated with reduced enzyme efficiency or DHEAS production in adults. Additionally, for all 29 positions, the modern human major allele is ancestral, and there is no evidence for a suite of novel adrenarche-extending variants. We conclude that the genetic foundation for extended adrenarche is shared between Homo sapiens and Neanderthals, and may have deeper hominin roots. Any phenotypic variation in adrenarche between modern humans and Neanderthals is more likely attributable to differential gene expression than to divergence in protein-coding sequences.

Mammal tails have long been recognized for their diversity of morphological form and function, however, there remains a substantial gap between the motivation to understand and emulate the various performance functions of the tail and what is known about tail anatomy. In this study, we were motivated to discover the anatomical foundations of the fast, whipping motions of the tail of the lesser Egyptian jerboa (Jaculus jaculus), which may aid in the quick changes of direction as the animal escapes from predators using ricochetal bipedal hopping. We employed microCT scans, dissections, and museum data to describe the musculoskeletal anatomy of the jerboa in comparison with the laboratory mouse (Mus musculus) and rat (Rattus norvegicus). While many aspects of tail anatomy are conserved across these species, the jerboa does possess unique characteristics such as an extremely long tail arising from caudal vertebral elongation, development of extensive dorsal musculature differentiated into lateral and medial components to increase points of skeletal attachment, and a novel anatomical feature - the bi-lobed cranial transverse process - which serves as a supernumerary dorsal tendon attachment site and possible brace to protect the ventral tendons and intrinsic muscles for a section of caudal vertebrae which likely experiences high mechanical stress.

Despite a common genome, males and females are remarkably different across species. Intralocus Sexual Conflict (IASC) arises when selection favours different trait values in the sexes but a shared genetic architecture constrains their evolutionary divergence. IASC is typically studied at the single-trait level, and genome-wide architectures can obscure conflict localized to specific genomic regions. We investigated IASC in the multivariate and local genetic architecture of 17 human metabolic traits and lifetime reproductive success by estimating the sex-stratified additive genetic (co)variance matrix (Gmf) and local cross-sex-cross-trait genetic correlations. Genome-wide, between-sex covariance matrix (B) showed only 1.15% asymmetry, and sexually concordant (SC) genetic variation exceeded sexually antagonistic (SA) variation by 18.5-fold. Indirect evolutionary response to SC selection was 1.6-fold stronger than direct responses to SA selection. In contrast, local analyses revealed a more heterogeneous structure - among the regions with significant local correlations, 12% were consistently SA, showing opposite fitness effects in the sexes, with functional enrichment in WNT signalling, while 20% exhibited both SA and SC effects . Our results demonstrate that in human metabolic traits, SC variation predominates genome-wide, whereas SA variation is concentrated in specific local regions, highlighting the importance of integrating multivariate and local frameworks in understanding conflict.

Why male primates invest in costly behaviours producing public goods remains debated, with two leading explanations, paternal care and reputation-based partner choice (RBPC). Using long-term data from four groups of wild vervet monkeys, we tested: (1) whether males show a bias in four protective "male services" (predator alarm calling, participation in between-group conflicts, leading river crossings and sentinelling); (2) which males contribute most; and (3) whether service provision predicts mating success during the mating season. We confirmed a male bias in all services. Consistent with the paternal care hypothesis, contributions were positively associated with past mating success, independently of rank, although potential fathers did not contribute more than non-fathers. Among non-fathers, service provision varied with rank, suggesting that newly immigrated males adjust their behaviour according to competitive state. Crucially, variation in alarm calling and between-group conflicts predicted future mating success, with between-group conflict emerging as the strongest and most consistent predictor of mating success across years and within mating seasons, whereas rank, tenure and social integration added little explanatory power. In contrast, sentinelling and leading river crossings did not reliably translate into mating benefits. Our findings indicate that male services are shaped by multiple selective pressures operating across different male career stages and that some forms of public goods provision function as signals of quality and cooperativeness to females. By directly linking cooperative investment to mating outcomes in a wild primate, this study provides rare empirical support for reputation-based partner choice beyond humans and highlights female choice as a potentially important force in the evolution of cooperation.

BackgroundGray is one of the relatively rare coat colors in donkeys. The Hetian Gray donkey is a distinctive indigenous breed from the Xinjiang Uygur Autonomous Region of Northwestern China, characterized by progressive hair depigmentation with aging while retaining dark skin pigmentation. However, the genetic basis underlying this unique gray coat color phenotype remains unclear. ResultsTo elucidate the genetic basis, we conducted whole-genome resequencing of Gray and non-Gray donkeys. Genome-wide selection signature analyses identified a candidate region on chromosome 15. Subsequent fine-mapping using mass spectrometry-based genotyping of 42 loci refined the candidate interval and revealed a SNP within intron 2 of the ASIP gene, located in a genomic fragment with highly similar sequences, showing complete association with the gray coat color. Association analysis in an expanded population further confirmed a strong correlation between this variant and the gray phenotype. Gene expression analyses also supported the role of ASIP in regulating pigmentation in donkeys. ConclusionsThese findings identify a genetic determinant of gray coat color in donkeys and provide new insights into the molecular mechanisms underlying age-related depigmentation in domestic animals.

Late Pleistocene sea-level fluctuations and glacial corridors intermittently connected the island of Ireland to Britain and continental Eurasia, shaping patterns of megafaunal dispersal and occupation. However, the scarcity of genetic data from ancient Irish fauna has limited our understanding of their demographic histories and relationships to continental populations. To address this, we generated sedimentary ancient DNA sequences from Castlepook Cave in southwest Ireland, detecting twelve ancient taxa, including two without zooarchaeological records at the site. We recovered the first mitochondrial sequences from Irish cave hyenas and woolly mammoths, providing new insights into their maternal population history: cave hyenas carried mtDNA haplogroup A1, previously identified in Late Pleistocene European populations, while woolly mammoths belonged to clade III/B2, which was replaced in Europe at approximately the same time. We also identify two mitochondrial clades (1b and 2), consistent with pre-Last Glacial Maximum mitochondrial lineage turnover in Ireland. Together, these results indicate that Irish megafauna formed a biogeographical continuum with their conspecifics in Europe, extend the known geographical range of several mtDNA lineages, and support faunal connectivity at the northwestern edge of Europe prior to the Last Glacial Maximum.

Remote island systems with small landmasses and reliable estimates of human occupancy are ideal model systems to disentangle the roles of global climatic changes and local human occupation on biota. Here, we used mitochondrial and nuclear genomic data from five endemic Fijian Lasioglossum (Homalictus) bee species to infer changes in effective population size (Ne). These ground-nesting bees are native, with non-specialised floral visitation habits, and distributed across the elevational gradient. All lowland species and populations showed strong signals of increasing Ne that correspond to the timing of human occupation of Fiji, but not Holocene climatic change. Highland populations, with greater isolation and present in regions less affected by anthropogenic impacts, did not show evidence of recent rapid increases in Ne. Population expansion rates across the elevational gradient differed between taxa, with significantly earlier and larger increases in predominantly lowland species than those with more restricted ranges in the highlands. This is consistent with the movement of people inland from coastal regions and into montane elevations of the island, and corresponding landscape changes that benefit the ecology of these bees. Specific life history traits of these bees, combined with substantive clearing of forest cover and floristic changes at lower elevations, has likely increased nesting opportunities and abundance of invasive floral resources. Our findings contrast with recent evidence that human occupation of Fiji has resulted in decreased ant biodiversity and raise the paradoxical possibility that human-mediated environmental changes may benefit some native montane tropical insect faunas.

Fossil data is crucial to construct phylogenetic time trees, which serve as the basis to test a wide range of evolutionary hypotheses. While the fossil record is known to be incomplete, modern stratigraphy provides predictions of the structure of the fossil record as expressed by gap location and duration. Advances in phylogenetic model development allow us to propagate this information into Bayesian phylogenetic inference in the form of priors on time-variable fossil sampling. However, the impact and role of stratigraphic architectures on time tree inference has so far remained unexplored. We introduce a novel simulation framework that combines realistic stratigraphic forward models with phylogenetic simulations. Using this framework, we examine (1) how stratigraphically plausible model violations of fossil sampling due to gaps affect total-evidence inference under the fossilized birth-death model and (2) if stratigraphic knowledge on gap duration and timing improves inference when incorporated in priors on fossil sampling. We find that total-evidence analysis is robust to stratigraphically plausible distribution of gaps in disparate stratigraphic architectures, with results being instead dominated by the number of morphological characters. Surprisingly, incorporating information on prominent gaps in the stratigraphic record does not improve phylogenetic inference. Our results suggest that phylogenetic inference is robust to model violations introduced by stratigraphic gaps over short timescales, with results being dominated by a priori known data availability constraints such as morphological character matrix size. This research establishes the foundations for joint modeling of phylogenetic and stratigraphic processes and narrows the knowledge gap between paleontology, stratigraphy, and neontology.

The timing of reproductive processes can influence how mating interactions shape reproductive strategies. We examined how female age correlates with receptivity and preference in the jumping spider Habronattus formosus, a system characterized by elaborate male courtship and strong female choice. To test how receptivity changes across the post-maturation period, we paired females of different ages with males and quantified male courtship displays and mating outcomes. Females were not receptive immediately after maturation and instead exhibited receptivity primarily 2-3 weeks after maturation. We further found that the full male courtship predicted mating success in older females, whereas only the second stage of the courtship predicted success in younger females. The delay in female receptivity is consistent with expectations for strong female choice systems, where females may have a greater opportunity to evaluate male courtship before mating. These results highlight the importance of age-dependent changes in female reproductive state and suggest that the timing of receptivity may shape how different components of male courtship influence mating success.

Ancient environmental DNA (aeDNA) from permafrost, lake, cave, and marine sediments provides a rich source of genetic data that captures broad perspectives of past biodiversity. Accurate dating is crucial for discovering ecologically relevant patterns from aeDNA, and molecular clock dating would allow for sample ages to be estimated from the recovered genetic material itself instead of the geological components. However, the fragmented and damaged nature of short-read ancient DNA (aDNA) from multiple taxonomic sources poses significant challenges and has limited this dating approach for aeDNA. Here we developed ratePlacer, a phylogeny-based method for analyzing aeDNA that can combine information from many short reads in a sample while accounting for DNA damage to provide maximum likelihood estimates of sample ages. Simulations demonstrate that ratePlacer accurately dates samples even under the fragmented, damaged conditions characteristic of aeDNA and outperforms Bayesian tip-dating approaches for taxonomically mixed samples commonly found in aeDNA. Yet age estimates from re-dating Kap Kobenhavn varied across taxa, highlighting the difficulty of molecular clock dating in aeDNA. This dating also revealed elevated G[->]T and C[->]A mismatches consistent with oxidative damage. These patterns reveal aDNA damage beyond deamination and that remains understudied, suggesting that aeDNA should be carefully evaluated in genomic and evolutionary analyses. The new dating method, ratePlacer, extends molecular clock dating of aDNA from single-specimen to pooled environmental DNA data, where traditional methods struggle.

Social hierarchies are a nearly universal feature of animal groups, but whether dominance reflects a single generalized trait or a collection of context-specific competitive abilities remains unclear. Here, we assess social hierarchy in three strains of laboratory mice (BALB/c, C57BL/6, and Shank3B knockouts) of both sexes using three established paradigms: the tube test, the warm spot assay, and the void spot assay. Hierarchies emerged in all strains and both sexes across all three assays, but how animals established rank differed markedly by strain and sex. In the tube test, Shank3b-/- knockout females, but not males, lacked the winner effects seen in wild-type mice, indicating that the ability to build a winning streak depends on social recognition in a sex-specific manner. In the warm spot assay, females formed stronger hierarchies than males, particularly among mice on a C57BL/6 background, with high-ranking females actively displacing others from the warm platform. In the void spot assay, BALB/c mice of both sexes frequently displayed territory-marking behavior, a pattern that was less common in the other strains. Overall, individual rank rarely generalized across domains, despite high trial-to-trial repeatability for individuals within each assay. Together, these findings indicate that mice behave as dominance specialists rather than generalists, with strain- and sex-specific strategies for establishing rank in different competitive contexts, suggesting that distinct neural circuits likely underlie these separable components of competitive ability.

Pathogen-mediated disease burden imposes some of the strongest selective pressures on human populations, shaping genetic variation in immune-related genomic regions. Because historical, cultural, ecological, and environmental factors can influence disease burden, associations between these factors and immune-related genomic signatures can provide insight into how pathogen-mediated selection varied among populations. To investigate these associations and identify factors potentially contributing to historical disease burden, we analyzed worldwide variation in the Major Histocompatibility Complex (MHC). We evaluated the relationships between global patterns of genetic variation and a comprehensive set of historical-cultural, ecological, and climatic variables compiled in DisECCO (Disease, Ecological, Cultural, and Climatic Origins dataset). This dataset includes factors potentially associated with historical disease burden, including the timing of major cultural transitions, historical climatic conditions, exposure to domesticated animals, and pre-industrial pathogen stress. Contrary to genome-wide expectations, MHC diversity and variation were not well explained by geographic distance from Africa, suggesting that neutral demographic processes play a limited role in shaping MHC variation. Instead, MHC genetic variation was significantly associated with climatic variables, domestication exposure, and the time lag between the onset of the Neolithic and urbanization, with the strongest explanatory models identified for MHC Class II. These results indicate that cultural transitions and environmental conditions have played a central role in shaping immune-related genetic variation, and may have contributed substantially to changes in disease burden. Overall, our findings show that integrating ecological and cultural-historical variables with genomic data can help explain global patterns of genetic variation in humans.

Recent evidence of extensive call sequence use in non-human primates has led to the theory that syntax evolved to mitigate the constraints of their genetically fixed repertoires, before vocal production learning later emerged in humans. However, evidence of similarly extensive sequence repertoires in an open-ended vocal production learner--the Western Australian magpie (Gymnorhina tibicen dorsalis)--offers a unique opportunity to explore potential alternative pathways to syntactic communication. Our previous work revealed fledgling magpies learn group-specific repertoires of structured call sequences from their social contacts, with more sociable individuals acquiring larger repertoires earlier in development. Notably however, the individual vocal segments that combine to form their calls and call sequences were shared across groups and emerged as early as the first week post-fledging--suggesting the underlying vocal elements may not be learned. Here we utilised acoustic neighbourhood-based dimensionality reduction to compare clustering patterns of vocal segments across magpie fledgling developmental stages, and between fledglings and adults. We found no evidence of acoustic development over time, and no significant distinction between fledgling and adult productions of the same vocal segments. The same coarticulatory effects--where a vocal element is produced differently when combined with another--and geographic variation established previously in adults were supported in fledglings too. These findings support that the vocal building blocks underpinning magpie call sequences are innate, suggesting usage learning better explains how fledglings learn to combine calls. In a species capable of open-ended production learning, this suggests learning to combine existing signals may be more adaptive than productively learning new ones. Rather than evolving solely to compensate for genetically fixed repertoires, syntax may have evolved as a flexible, convergent solution to the various challenges of expanding communicative capacity--whether due to genetic constraints, cognitive limitations or the cost of establishing new meaning in novel signals.

The study of historical documents from the original excavation led by J. M. Cruxent, of museum collections, as well as results of fieldwork, are combined to provide a revised assessment of the archeologically site of Cucuruchu, which together with those of Taima-Taima and Muaco, document megafaunal lithic remains associated from the Late Pleistocene in northern Venezuela. Six vertebrate taxa are reported for Cucuruchu: five megaherbivores (a giant sloth, Eremotherium laurillardi; a glyptodontid, Glyptotherium cf. G. cylindricum; a macrauchenid, cf. Xenorhinotherium bahiense; a proboscidean, Notiomastodon platensis; and an equid, Equus sp.), and a turtle (Kinosternon sp.). Two fragments of El Jobo projectiles were found in situ lying among megaherbivores bones, while another unidentified projectile fragment was also found in the same layer. Although no direct evidence of human exploitation on these megaherbivores remains has been found, the presence of El Jobo projectiles in the undisturbed fossil-bearing layer some of which are associated with fossil remains, could be a suggestive of such interaction. Two new radiocarbon dates indicate ages of at least ca. 15.3-16.0 kyr cal BP for the fossil-bearing layer, contrasting with the mid to late Holocene age originally reported. Cucuruchu is a window to the past and supports the presence of hunter-gatherers of the El Jobo culture in the region at the end of the Pleistocene.

Sex chromosomes are central to many ecological and evolutionary processes. Evidence has accumulated that sex chromosome systems vary extensively in age, turnover and transitions, motivating renewed efforts to study the diversity of sex chromosome systems across the tree of life. However, successful genomic detection of sex chromosomes depends on several factors, including the size and divergence time, background genetic diversity, and the number of sequenced females and males. In addition, technical challenges associated with sequencing and analysing the sex-limited Y/W chromosome remain. Here, we present PhaseWY, an automated Snakemake pipeline that uses whole-genome sequencing data from multiple female and male individuals to identify sex-chromosomal regions and extract the corresponding Y/W sequences. PhaseWY (i) detects sex differences in alignment depth, (ii) applies read-based and statistical haplotype phasing, (iii) identifies sex-linked regions using haplotype clustering, and (iv) subsets autosomal, X/Z- and Y/W-linked variants for downstream analyses. We applied PhaseWY to simulated data to benchmark factors influencing sex-linkage detection and successful extraction of Y/W-linked variants. To demonstrate its practical utility, we further applied PhaseWY to the neo-sex chromosome system in Alauda larks (Alaudidae) and performed a range of downstream analyses demonstrating the scope of applications of the PhaseWY output. We conclude that PhaseWY provides an easy-to-use and reproducible tool for population-genomic analyses in non-model organisms, with particular importance for advancing our understanding of sex-chromosome evolution.

Recent studies have characterised the mutational profile across multiple mammalian species, highlighting substantial differences across lineages. However, none of these studies investigated whether mutation profiles and geography are significantly correlated. In this study, we present a multi-genome alignment spanning several Capra taxa, reconstruct the ancestral genome of Capra hircus and use it to characterize the mutational profiles across multiple Capra species by using the 1000 genomes VarGoats dataset. Results confirmed that the scale of differences among Capra species largely reflects their phylogenetic relationships, in particular with the Bezoar being genetically closer to domestic goats than to other wild species. Subsequently, we correlated the mutational profile and the geographical origin of the different individuals. In particular, ACG>ATG changes have the strongest correlation with longitude (r = -0.79, P-value = 3.02*10-204), while TCA>TGA are strongly correlated with latitude (r = -0.51, P-value = 4.30*10-63). We highlight how sequential dinucleotide mutations (SDMs) place cosmopolitan breeds closer to the sampling location, rather than the country of origin, showing how the recent relocation of cosmopolitan breeds to new continents is reshaping the genome of these animals. Finally, we used the mutational profile to predict the coordinate of origin of each animal in the dataset. In conclusion, we show the important role that geography had in shaping the genomes of domestic goats.

The second plague pandemic (early 14th-early 19th centuries), which was caused by Yersinia pestis, had a profound demographic, socio-economic and cultural impact across Eurasia and North Africa. Many regions in Europe and the Middle East are estimated to have lost 40-60% of their human populations, with some areas suffering even higher mortality. Whether exposure to Y. pestis drove strong positive selection on protective genetic variants in the human genome, and how it shaped migration patterns, remains debated, despite several recent studies based on ancient DNA. Here, we analyse a markedly larger, higher coverage, and geographically diverse dataset based on shotgun sequencing of genomes from 529 ancient individuals to a mean depth 8.8x dating to either before or after the arrival of the pandemic at three sites in northern Europe: Trondheim (Norway), Lund (Sweden) and Vilnius (Lithuania). Genome-wide scans for signatures of selection provide no evidence for strong positive selection acting on specific genetic variants driven by Y. pestis exposure: we neither replicate selection signatures reported by previous studies nor identify new genome-wide significant candidates. However, for all three sites, we observe evidence for a reduction in long-range immigration, indicated by a drop in the diversity of ancestry that followed the arrival of Y. pestis and broadly coincided with the end of the Viking Age, Christianisation and the onset of the Little Ice Age. Our results shed important light on the demographic impact of major sociohistorical changes that occurred during the late Medieval period in Scandinavia and the Baltic region and link Christianisation to increased diversity in ancestry before the pandemic.